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We investigate the growth of connectivity in a network. In our model, starting with a set of 
disjoint nodes, links are added sequentially. Each link connects two nodes, and the connection rate 
governing this random process is proportional to the degrees of the two nodes. Interestingly, this 
network exhibits two abrupt transitions, both occurring at finite times. The first is a percolation 
transition in which a giant component, containing a finite fraction of all nodes, is born. The second 
is a condensation transition in which the entire system condenses into a single, fully connected, com- 
ponent. We derive the size distribution of connected components as well as the degree distribution, 
which is purely exponential throughout the evolution. Furthermore, we present a criterion for the 
emergence of sudden condensation for general homogeneous connection rates. 
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Networks are sets of nodes that are connected by links. 
Models for the evolution of complex networks fall into 
two general classes: evolving graphs where the number of 
nodes is fixed but the number of links grows, and grow- 
ing networks where the number of nodes and the number 
of links both grow [H-Q- The classic evolving random 
graph model, where pairs of randomly selected nodes are 
repeatedly connected by links, captures the nucleation of 
a giant connected component with a macroscopic num- 
ber of nodes J,, J3] . The preferential attachment model 
of network growth, where newly added nodes connect 
to existing nodes with probability that is proportional 
to the degree, yields the broad degree distributions with 
power-law tails that characterize many real-world com- 
plex networks [HI ■ 

In this letter, we study random graphs that evolve 
according to a preferential attachment mechanism. In 
our model, the number of nodes is fixed but the num- 
ber of connections grows. The "popularity" of each 
node, as measured by the degree, governs the connection 
process. This fusion between the two seminal network 
models of random graphs and preferential attachment, 
is inspired by Facebook, the immense online network 
of cyber- friends Q. In Facebook, new connections are 
formed via friendship requests from one member of the 
network to another; upon acceptance, the two become 
friends. Preferential attachment implies that members 
seek and accept friends according to popularity. 

In our model for the popularity-driven growth of con- 
nectivity in a network, the system consists initially of a 
set of disjoint nodes. We consider the natural situation 
where the degree of a node, defined as the number of ex- 
isting connections, controls the rate by which two nodes 
connect with each other. Specifically, a node with degree 
i and a node with degree j connect at a rate that is lin- 
ear in i and linear in j as well. We find that the degree 
distribution remains purely exponential throughout the 
evolution. Thus, the preferential attachment mechanism 
leads to opposite results in fixed and in growing networks: 
the tail of the degree distribution is narrow in the former 




FIG. 1: Connection and subsequent aggregation. The new 
link (dashed line) connects two nodes with degrees i = 3 and 
j = 2. As a consequence, two connected components with 
sizes I — 4 and m = 3 merge. 



but broad in the latter. 

Our main result is that the system undergoes two 
finite-time transitions. In the first transition, occurring 
at time tg, a giant connected component nucleates. The 
giant component contains a finite fraction of all nodes. In 
the second condensation transition, the giant component 
takes over the entire system. In contrast with classical 
random graphs, the entire network condenses into a sin- 
gle component at a finite time tc- We obtain analytically 
the degree distribution and the size distribution of con- 
nected components when the rate of connection between 
nodes with degree i and nodes with degree j equals the 
product {i + + 1). From these two distributions, we 
obtain the finite transition times: 



tg = 1/3 



and 



tr. 



1. 



(1) 



The rest of this letter includes two parts. In the first, we 
study how the degree of a node grows as a result of the 
networking process, and in the second, we investigate 
how the size of a connected component increases as a 
result of the very same process. 

In our evolving random graph model, the network con- 
sists of a set of disjoint nodes at time t = 0. The net- 
work becomes connected through the sequential addition 
of links: a node with degree i and a node with degree j 
connect with rate Cij. As a result, the degree undergoes 
the augmentation process (figure [T|) 



+ 1,3 + 1). 
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Since the two nodes are interchangeable, the connection 



rate is symmetric, C^^ 



We note that the clas- 



sic random graph corresponds to the uniform connection 
rate dj = 1 [i,|9l|. 

Throughout this study, we implicitly take the infinite 
system size limit. Let nj{t) be the degree distribution, 
that is, the fraction of nodes with degree j at time t. This 
quantity is normalized, J2j '^i = 1; ^^'^ it obeys the rate 
equation 



drij 
~dt 



(3) 



The initial condition is rij(O) = Sj^. The quantity Vj 
equals the total connection rate of nodes with degree j. 

By summing the evolution equations ([3]), we can ver- 
ify that the normalization is preserved, (d/dt) rij = 0. 
Similarly, by multiplying ^ with j and summing over all 
J, we find that the average degree, (j) = '^jjnj, grows 
according to 



(4) 



This equation reflects the connection process (HJ. Since 
every link connects two nodes, the average degree conve- 
niently yields the total density of links, L, via the simple 
relation 2L = (j). 

In the preferential attachment model of network 
growth, the attachment rate is linear in the degrees of 
the existing nodes. In this "rich-get-richer" mechanism, 
the connection probability is proportional to the popu- 
larity. Hence, to model popularity-driven networking in 
evolving graphs, we restrict our attention to rates that 
are linear in both i and j, Ci,j ~ (i+a)(j + a). Often used 
in preferential attachment, the offset a allows us to start 
with a set of disjoint nodes, the natural initial condition. 
Here, we investigate the rate 



(5) 



We verified that the qualitative behavior, including the 
the sudden condensation transition, extends to all a . 

By substituting the rate (O into equation dU, we see 
that the average degree obeys the closed equation 



dt 



= {{j)+iy 



(6) 



We obtain the average degree by solving this equation, 
subject to the initial condition 0)|j_q = 0, 



1-t' 



(7) 



when t < 1. Hence, the average degree diverges in finite 
time. 

From the average degree ([7]), we find that the connec- 
tion rate Vj is linear in the degree, vj = (j -I- 1)(1 — t)~^. 



Therefore, the degree distribution obeys the linear evo- 
lution equation 



(1 



, drii 
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' dt 



(8) 



We solve these equations recursively, starting with the 
initial condition nj{0) = Sj^o, and obtain uq = 1 — t, 
ni = (1 — t)t, n2 — {1 — t)P, etc. Generally one finds 
that the degree distribution is purely exponential. 



(l-t)t^-. 



(9) 



Surprisingly, the degree distribution vanishes in finite 
time: Ujitc) = for all j with tc — 1. Consequently, 
the degree of a node is finite if and only ii t < tc- This 
behavior is a signature of the continuous condensation 
transition that occurs at the condensation time tc- As 
shown below, the entire network condenses into a single, 
fully connected, component at this time. 

The difference between growing and fixed networks 
is remarkable. Networks that grow by preferential at- 
tachment have broad degree distributions with power-law 
tails [13, im E^-O , but in a fixed network, preferential 
attachment leads to a narrow degree distribution with 
exponential tail. On its own, the rich-get-richer mecha- 
nism does not generate broad tails, but rather, it is the 
combination of a growing network and preferential at- 
tachment that leads to a broad distribution of degrees. 
We note that exponential degree distributions occur in a 
variety of complex networks including power, transporta- 
tion, and social communication networks [iq . [iTj . 

Our evolving graph consists of multiple connected com- 
ponents ("clusters") which undergo binary aggregation 
as a result of the connection process Symbolically, 
such an aggregation process [H, [l^ can be represented 
as (figure [T|) 
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(10) 



where / and m are the number of nodes in the two merg- 
ing clusters. The rate of aggregation between a finite 
cluster of size I and a finite cluster of size m is 



if/,™ - (3;-2)(3m-2). 



(11) 



To obtain this rate, we note that a cluster with k nodes 
has k — I links. (This relation is valid for clusters with 
tree structure, and indeed, nearly all finite clusters are 
trees in the infinite size limit 7].) Equation (jlip follows 
from the connection rate ([5]) , together with the fact that 
the sum of the degrees in the cluster equals twice the 
number of links. 

The total density of clusters, c, follows from the total 
density of links, L, when all clusters are finite. Since 
every link reduces the number of clusters by one, we have 
c(t) = l~L{t). Using 2L = (j), we conclude c = l-(j)/2, 
and from ([7]), we obtain the cluster density 



c(t) 



2-3i 
2(1 - t) 



(12) 
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As shown below, this relation holds when t < tg with tg 
given in ([T]). 

The density Ck {t) of clusters with size k at time t obeys 
the master equation 

% = 5 II (3;-2)(3m-2)QC,„ - (j + l)(3fc-2)cfc, (13) 



l-\-m—k 



where (j + 1) = (1 — t)^^. The initial condition is 
Cfe(O) — Sk,i- The gain term is nonlinear, and it directly 
reflects the aggregation process pH)) with the product 
rate pT|) . The loss term, on the other hand, is linear in 
the cluster-size density. The loss rate equals the prod- 
uct between the sum of all degrees (shifted by one) in 
the cluster, and the average degree (shifted by one) of all 
nodes in the entire system. In this form, the master equa- 
tion is valid at all times, whether all clusters are finite or 
whether macroscopic clusters exist as well (20l |. 

Let M„ — '^"cfc be the nth order moment of the 
size distribution. The zeroth moment gives the total clus- 
ter density, Mq = c, and the first moment. Mi, yields the 
total mass of finite components. Furthermore, the second 
moment satisfies 



dt 



= (3M2 - 2Mi)2 

+ (3M3 - 2A/2) [(3A/i - 2Mo) - (j + 1)] , 



subject to the initial condition M2{0) = 1. We ob- 
tain this equation by multiplying (|13p by and sum- 
ming over all k. Let's assume that finite clusters contain 
all the mass. Mi = 1. In this case, the cluster den- 
sity (IT2|) gives 3Mi — 2Mo = (1 — i)~"^, and as a conse- 
quence, the second moment satisfies the closed equation 
dM2/dt = {MI2 - 2)2. Therefore, 



l-2t 
1 - 3t 



(14) 



for t < 1/3. The divergence of the second moment 
shows that a percolation transition [2l[ occurs at time 
tg = 1/3 as stated in ([1]). The critical density of links 
Lg — L{tg) — 1/4 is smaller than the value Lg = 1/2 cor- 
responding to the classic random graph ^] . 

When t < tg, the system is in a non-percolating phase 
and finite clusters contain all the mass. Mi = 1. Other- 
wise, the system is in a percolating phase, where a macro- 
scopic cluster, the giant component [2^ . contains a finite 
fraction of all nodes and Mi < 1. 

To find the mass of the giant component, we analyze 
the cluster-size distribution. First, we solve recur- 
sively for small clusters. We can confirm that the den- 
sity of minimal-size clusters ("monomers") equals the 
fraction of isolated nodes, ci = riQ = 1 — t. The densi- 
ties of "dimers" and "trimers" are C2 = ^t{l — t)^ and 

= t^{l ~ t)^. We therefore expect the general form 



Ck 



Ak t''-\l-t) 



2fc-l 



(15) 



By substituting this expression into (|13p . we obtain 
recursion equation for the coefficients Ak 



2{k-l)Ak= ^ {il-2){im-2)Ai 



-lA,m 



(16) 



^va—k 



for fc > 1, subject to Ai = 1. We now make the trans- 
formation Bk = (3A; — 2)Ak- The coefficients Bk satisfy 
a second recursion relation 

2(fc-l)Sfe = (3fc-2) 

l+m—k 

for fc > 1, subject to -Bi = 1. From this recursion, we 
conclude that the generating function B{x) = J^k x'', 
obeys the differential equation = i_^g ■ We now 
integrate this differential equation subject to the con- 
straint B = 1, and find that B{x) satisfies 
the cubic equation 



B(l -B)= 



(17) 



We obtain the coefficients Bk by using the Lagrange 
inversion method, 



Bk — 



1 

2Td 

1 
27ri 

1 
2Td 

1 

2Td 



B 



dx 



B[l - B){1 - 3B) 

_ B)2(fe+1) 

1 - 3B 



dB 



3fc- 1 
fc - 1 



B^{1 - Bf^^ 

2k + n 
n 

3k - 2 
fc-2 



.n=0 



dB 



dB 



(3fc-2)! 
fc!(2fc- 1)!' 



In the second line, we used equation (jl7p . and replaced 
the integration over x with integration over B, by using 
dx = {I ~ B){1 — 3B)dB. Finally, we utilize the identity 
(1 - S)-™ = X;r=o to determine the coeffi- 

cients Ak, 



Ak = 



(3fc-3)! 



(18) 



fc!(2fc- 1)!' 

{2,1,2,6,22,91,...} also 



The integer sequence 2Ak 
arises in planar maps [2^ [2^ . 

The tail of the cluster-size distribution is a product of 
a power law and an exponential. 



for fc 



1 



00, with k^ ^ 



, ^-5/2 g-fe/fc. 



In 



tg(l-ta)" 



(19) 



This tail fol- 



lows from the Stirling formula and the coefficients (|18p . 
Hence, at the percolation transition, the distribution is 
purely algebraic, Ck{tg) ~ fc"^/^ d^]. Otherwise, the tail 
is power law at small scales, k <^ k^., but exponential at 
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FIG. 2: The mass of the giant component, g, and the total 
density of connected components, c, versus time t. 



large scales, k k^,. The characteristic scale diverges, 



k^ — '2j{,^ ^p) ' ^ 

Note that the time dependence of the cluster-size den- 
sity ([T5|) enters primarily through the quantity t{l — t)"^ 
as Ck oc [t(l — i)^]'^. This observation tells us that for an 
auxiliary time r, that is connected to the physical time t 
via the cubic equation 



r(l-r)2=i(l-^f 



(20) 



the densities Ck{t) and Cfe(T) are related by the duality 
relation 



1-t 

1-T 



(21) 



When t < tg, equation ([^0)) has only the trivial solu- 
tion T — t, while for tg < t < tc there is additionally 
a second nontrivial root t < tg. Consequently, for all 
tg < t < tc we can choose the nontrivial root t < tg 
of (PO)) and then, the duality relation (PT|) conveniently 
specifies the cluster-size distribution in the percolating 
phase in terms of its counterpart in the non-percolating 
phase. By multiplying (1211) by k and summing over all 
k we see that Mi(t)/(1 — t) = Mi(t)/(1 — r), and since 
Mi(t) = 1, we immediately get the nontrivial first mo- 
ment in the percolating phase, Mi = (1 — t)/(l — r). In 
summary, the mass of finite clusters is 



1 



< t < i 



9' 



t> tr. 



(22) 



The mass of the giant component, g, equals the comple- 
mentary mass, g = 1 — Mi (see Fig.[2]). Post-percolation, 
this quantity grows linearly, g{t) ~ 3(t — tg), as t tg, 
and furthermore, 1 — g{t) :^ tc — t a.s t 'f tc- At time 
tc — 1, the giant component takes over the entire graph, 
and hence, the system condenses into a single compo- 
nent. Remarkably, the cluster-size density and the frac- 
tion of nodes with finite degrees vanish simultaneously, 
Ck{tc) = rijitc) = 0. 



For completeness, we mention that the generating 
function, C{z,t) — '^k{t) , can be written as an ex- 
plicit function of time using the hypergeometric function. 



Q{z,t) = 



1 



3t(l -t) 



F 



- 1 



Moments of the distribution can also be written explicitly 
[11]. The density c{t) = Q[z = l,t) is plotted in figure 
m and one can verify that this quantity matches ((T2|) in 
the non-percolating phase. 

An unusual feature of our network connection process 
is that both the percolation transition and the conden- 
sation transition occur at finite times. Let us consider 
the uniform connection rate = 9, corresponding to 
the classic random graph. For this evolving graph, the 
percolation time is finite, tg — 1/9, but the condensa- 
tion time is divergent, tc — oo (in a finite system, the 
condensation time is logarithmic in the total number of 
nodes) 01 • Interestingly, the corresponding aggregation 
rate, Ki^m = (3Z)(3m), is Zander than the aggregation rate 
pTjl . yet the latter, smaller, rate produces faster conden- 
sation! Of course, the aggregation rate Ki^rn corresponds 
only to finite clusters, and it does not apply to the gi- 
ant component. The rich- get-richer mechanism acceler- 
ates the rate by which the giant component engulfs finite 
components, and this feature is ultimately responsible for 
the finite-time condensation. 

We also considered a general class of evolving graphs 
with homogeneous connection rates, and we now briefly 
discuss the purely algebraic rates 



(U) 



(23) 



The case a = corresponds to ordinary random graphs, 
and the case a = 1, to popularity-driven networking. 
From (ID), the average degree grows according to 



dt 



if 



{3? 



(24) 



where we assumed the scaling behavior ^ ■ 

From this rate equation, we find the scaling behaviors 



^l/(l-2a) 
„ const. X t 



a < 1/2, 

e-"'" a = 1/2, (25) 
(t,-t)-i/(2"-i) l/2<a<l. 



When a < 1/2, the average degree grows algebraically 
with time. Furthermore, the percolation time is finite, 
but the condensation time is infinite. In the marginal 
case a = 1/2, the degree grows exponentially with time. 
When 1/2 < a < 1, the percolation time is finite, and the 
condensation time is finite as well. The average degree 
diverges as the condensation transition is approached. 

Finally, when a > 1, the scaling assumption used in 
p4| is no longer validj and condensation becomes instan- 
taneous: tc = tg = 27]. This phenomenon is analogous 
to the instantaneous gelation that occurs in aggregation 
[28l [29| and exchange-driven growth ^] . 
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In conclusion, we studied the growth of connectivity in 
an evolving graph. In our model, the degree controls the 
connection process as the probability that a node forms 
a new connection is proportional to the total number of 
its existing connections. We find that the system under- 
goes two continuous transitions, both occurring at finite 
times. In the first percolation transition, a macroscopic 
connected component nucleates. In the second conden- 
sation transition, the network becomes fully connected. 
We obtained analytically the degree distribution and the 
size distribution of connected components. We find that 
in a fixed network, the degree distribution is exponential, 
and therefore, the average degree fully characterizes the 
entire distribution. 

Our theoretical analysis relies heavily on the fact that 
the connection rate is linear in the degree. For such con- 
nection rates, the rate by which components merge is a 
bilinear function of size. Consequently, the size distribu- 
tion of connected components is analytically tractable. 



Linear connection rates also have computational advan- 
tages. Popularity-driven networking can be studied using 
a convenient random link algorithm. In this numerical 
simulation, we choose two links at random. For each 
link we pick one of the two nodes connected to it, and 
then, connect these two nodes. With this implementa- 
tion, we select nodes with probability that is proportional 
to the degree 0], and this procedure can be generalized 
to arbitrary linear rates. The random link algorithm has 
optimal efficiency as the computational cost is linear in 
the number of links, and it can be used to study addi- 
tional properties of our evolving graph. Among a host 
of possibilities, we mention structural properties of con- 
nected components such as their cycle structure [12] , and 
extremal properties such as statistics of the largest con- 
nected component. 
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